Quantum size effect in Pb(lOO) films: the role of symmetry and implication for film 

growth 
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We show from density-functional calculations that Pb(100) thin films exhibit quantum size effect 
with a bilayer periodicity in film energies, film relaxations, and work functions, which originate from 
different symmetry of the stacking geometry of odd and even layer films. The bilayer periodicity of 
the film energy is argued to survive on a semiconductor substrate, which should allow the growth 
of "magically" thick even-layer Pb(fOO) films. Furthermore, it is found that the quantum well 
states in a simple metal film can be classified into a-bonded and 7r-bonded states, which quantize 
independently. 
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The so-called "electronic growth" model of metallic 
thin film on a semiconductor substrate has attracted 
much scientific and technological interest Q . This model 
provides conditions for the growth of atomically-flat thin 
films or nano-structures, which are of significant techno- 
logical interest. The underlying physical mechanism be- 
hind "electronic growth" is the quantum size effect (QSE) 
on the film energy introduced by quantum well (QW) 
states in the film. In addition to the film energy, the 
energy gain associated by the charge transfer occurring 
at the film semiconductor interface needs also to be con- 
sidered to determine the film stability. The total energy 
of the system may then result in "critical" thicknesses in 
which the film is stable when the number of atomic layers 
is above or below a critical number and "magical" film 
thicknesses in which the film has a pronounced stabil- 
ity for certain number of atomic layers. The "electronic 
growth" model has been very successful in understanding 
observed "critical" thicknesses of alkali metal and noble 
metal overlayers grown on a GaAs surface 0, and the 
observed "magical" thicknesses of Pb(lll) films or islands 
on various substrates 0, HUH. 

The nature of observed QW states and the associated 
QSE in metal film have so far been well understood from 
simple models based on bulk band structure and Fermi 
surface nesting^, 0, 0] In general, QW states are 
formed by the quantization of the valence electrons con- 
fined in the perpendicular direction of the thin film, as 
described by the phase accumulation model. The ener- 
gies of the QW states are then obtained from the energies 
of the valence bulk bands at the quantized wave vec- 
tors. This simple model has been widely used in analysis 
of photoemission measurements of Q W states 0, IS ■ 
The QSE on the film energy and other physical proper- 
ties arises from the consecutive occupation of QW states 
with increasing film thickness. An enhanced density of 
states at the Fermi surface by QW states is then obtained 
from extremal spanning (nesting) vectors q of the Fermi 
surface that are perpendicular to the film, resulting in 



QSE period Xqse = 2ir/q for the film energy. For a 
free-electron model of the bulk bands, corresponding to 
a spherical Fermi surface, q = 2fcp where fcp is the Fermi 
wave vector and Xqse = it /hp- 

This simple model for QSE has been supported by 
first-principles, electronic structure calculations for var- 
ious systems. Schulte ^3 found that the electron den- 
sities, potentials, and work functions of a thin jellium 
film exhibited QSE with a period tv/Icf in accordance 
with a free-electron model. Feibelman ^3 and Batra 
et. ai. 12] found a QSE in the surface energy and sur- 
face relaxation of both Al(lll) and Mg(OOOl) thin films. 
Wei and Chou[l3 found that the QSE period for surface 
energies and work functions of thin Pb(lll) films was 
consistent with a spanning vector in the Pb bulk band 
structure in [111] direction. The observed oscillations of 
film thickness and interlayer spacing had the same QSE 
period 0, H H El E3 So far the calculated and ob- 
served QSE periods can be simply understood from a 
nesting vector of the bulk band structure along the di- 
rection perpendicular to the film, which projects to the 
center of the film Brillouin zone (FBZ). 

In this letter we show from a density-functional the- 
ory (DFT) study of thin Pb(lOO) films that the symme- 
try of the stacking geometry of a thin film can have a 
profound influence on the QSE through QW states at 
the FBZ boundary. In particular, the film energy, relax- 
ations and work functions exhibit a bilayer periodicity 
that cannot be understood in terms of nesting vectors of 
the Fermi surface along the TX direction in the bulk BZ. 
This symmetry effect is a general phenomena for fcc(lOO) 
(and bcc(lOO)) films but should only be important for the 
film energy when the material also has a large band gap 
around Fermi level at the FBZ boundary. 

We have carried out DFT calculations of the total en- 
ergy, geometric and electronic structure of freestanding 
Pb(lOO) films using a pseudo-potential, plane-wave basis 
method within the local-density approximation. Techni- 
cal details can be found in Ref. |l7| . Spin-orbit coupling 
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FIG. 1: Film energy Eo per unit cell (left axis) and geo- 
metrical relaxations 3d and Sdi2 (right axis) of free-standing 
Pb(100) films with respect to the number of atomic layers. 
Eo of both the bulk truncated and the fully relaxed film are 
shown. Eo is obtained by subtracting the bulk contribution 
from the total energy of the film per unit cell. Sd and Sdi2 are 
the changes of the total film thickness and the first interlayer 
spacing from their truncated bulk values, respectively. Note 
that the energy and the relaxations exhibit a characteristic 
odd-even-layer alternation. 

was not considered. 

As shown in Fig. 1, the calculated film energies and re- 
laxations of the free-standing, thin Pb(100) films exhibit 
pronounced QSE. All these quantities exhibit a charac- 
teristic damped odd-even-layer alternation. The film en- 
ergies of the relaxed films have a much more pronounced 
bilayer oscillation than the bulk truncated film, showing 
that the QSE on the relaxations contribute to the oscil- 
lations of the film energy in a cooperative manner. Both 
the film thickness and the first interlayer spacing show 
very similar oscillatory behavior. Also the work func- 
tions of the relaxed films (not shown here) exhibit QSE 
and oscillate in an anti-phase manner relative to the os- 
cillations of the film energy (i.e, work function maxima 
appear at the film energy minima). 

A key issue is to understand why the QSE in thin 
Pb(100) films has a bilayer period. We first show that 
this periodicity cannot simply be understood in terms of 
nesting vectors of the Fermi surface along the TX direc- 
tion. The band structure around the Fermi level is dom- 
inated by bands derived from atomic p orbitals. Along 
the TX direction, corresponding to the perpendicular di- 
rection of the film, the Fermi surface has a nesting vector 
q = 1.367r/a corresponding to Xqse = 2.9<iioo where dioo 
is the inter-layer spacing. Thus a QSE period of about 
three (100) layers is expected in contrast to the calculated 
bilayer period for the QSE. 

As a first step to understand the calculated QSE pe- 
riod, we have studied the QW states of the isolated 2 
to 12 layer thick Pb(100) films. In Fig. 2 we show, as 
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FIG. 2: Band structure of a four-layer, free-standing Pb(100) 
film along high symmetry directions in the film Brillouin zone 
(FBZ). These directions are depicted in the irreducible part 
of the FBZ. The symmetry characters of the various states 
at the T point and along the symmetry lines in the FBZ are 
indicated. Note that the states along the Y line are two-fold 
degenerate in an even-layer film, and are nondegenerate in an 
odd-layer film. Fermi energy is set to eV. 

an illustrative example, the band structure of a bulk- 
truncated, four-layer Pb(100) film in the irreducible part 
of the film Brillouin zone (FBZ). The symmetry charac- 
ters of the bands are indicated. For each wave vector in 
the FBZ, the discrete set of bands corresponds to QW 
states. Note that the bulk s and p subbands are sepa- 
rated by a gap of about 4 eV. The number of QW (s-QW) 
states derived from the bulk s bands is given by the num- 
ber of layers (N) in the film, whereas there are 3N p-QW 
states derived from the bulk p bands. Only p-QW states 
cross the Fermi level and are relevant for the QSE. The 
p-QW states disperse either to higher energy or lower 
energy away from the T point depending on their orbital 
character being p z (cr)- or p^^-like. These states ex- 
hibit many "avoided crossings" accompanied with inter- 
change of their characters. A key observation about QW 
states of these films that cannot be understood from the 
bulk band structure is the lifting by odd-layer films of 
the two-fold degeneracy of QW states along the Y sym- 
metry axis in even-layer films. Note that the projected 
bulk band structure on this axis is two-fold degenerate. 

To understand the effect of the odd-even layer alter- 
nation of the degeneracy on the consecutive filling of the 
QW states, we show in Fig. 3 the calculated layer de- 
pendence of the energy spectra of p-QW states at high 
symmetry points in the FBZ. The p-QW states of ir bond- 
ing character at T and X points and p-QW states of a 
bonding character at the M point do not cross the Fermi 
level and are accordingly not shown 18]. The occupied 
bands at T can not accommodate all the p-electrons. In 
principle, it is necessary to consider the QW states at 
the whole FBZ, in particular at all the high symmetry 
points. 
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FIG. 3: Energy spectra of quantum well states derived from 
bulk p bands at the high symmetry points F, X and M in 
the film Brillouin zone as a function of the number of atomic 
layers in the film. At the V and X points only states with 
a bonding character are shown, whereas at the M point only 
states with n bonding character are shown. The states with 
different bonding characters at these k points do not cross 
the Fermi level and are not shown. Note that the splitting 
of the states around the Fermi level in panel (b) exhibit an 
odd-even-layer alternation. The Fermi energy is set to eV. 



Only the QW states at the X point have a bilayer peri- 
odicity in Fig. 3. It is clear that the QW states at the T 
point cross the Fermi level by a period of about 3 layers 
in good agreement with the perpendicular spanning vec- 
tor in the TX direction. Note that at the M point, the 
two- fold degeneracy of the 7r-like QW states are not lifted 
by odd-layer films. However, at the X point (and along 
Y symmetry axis) this degeneracy is lifted by odd-layer 
films corresponding to a bilayer periodicity. In partic- 
ular for states crossing the Fermi level, the splitting is 
anomaly large because of a large gap of about 1.5 eV in 
the projected bulk p band structure along the Y symme- 
try axis. The zone boundary band-gap in the thin film 
therefore exhibits an odd-even-layer oscillation, similar 
to that found in Si(100) thin films |Tg| . 

In general, QW states from different regions of the FBZ 
all contribute to the film total energy. In particular, the 
QSE period will be determined from the extremal per- 
pendicular spanning (nesting) vectors which determine 
the period of the enhanced DOS at Fermi-level. In case 



that QW states at zone center dominate, free-electron- 
like model is able to give reasonable prediction of the 
film stability, e.g., Pb(lll) film[3, Ag(100) film [2j. In 
the case of the Pb(100) thin film, the QW states at BZ 
boundary give rise to the periodicity of the film energy. 

The different degeneracies along the FBZ boundary in 
odd and even layer Pb(100) films is due to the symme- 
try of the layer stacking. Because the Pb(100) film (and 
the (100) film of any fee or bec metal) has an ABAB... 
stacking sequence, odd- and even-layer films have differ- 
ent symmetry. An odd-layer film has an mirror plane on 
the central layer. The point group is D^h, and the space 
group is symmorphic. For an even-layer film, there is no 
such direct mirror plane but a glide plane. The result- 
ing space-group D^ h is non-symmorphic. For the odd- 
layer film, the point group along Y axis is C^, and the 
bands of odd-layer film are always non-degenerate along 
Y. For the even-layer films, non-primitive translations are 
involved along this axis and a multiplier representation 
has to be considered ( see e.g., Ref. |21j). Only a single 
two-dimensional representation is found and the bands 
in even-layer films are all two-fold degenerate along the 
Y axis. 

So far the observed and calculated QW states of thin 
metal films has been successfully described by the bulk 
band structure and the phase accumulation model. How- 
ever, as illustrated by the results for the QW states of 
thin Pb(100) films, the energies of the QW states at the 
FBZ boundary of (100) film of any fee (or bec) metal are 
determined not only by the film thickness, but also by the 
symmetry of the layer stacking geometry. A mirror plane 
for odd-layer film and a glide plane for even-layer film 
result in non-degenerate and two-fold degenerate QW 
states along Y axis. This symmetry effect can not sim- 
ply be captured by the bulk band structure and the phase 
accumulation model. However, the size of the splitting 
across the Fermi level and the QSE in film energy and 
relaxation is material specific because it depends on the 
presence of a large band gap at the Fermi level at the 
boundary of FBZ. For instance, in the case of a Al(100) 
film, we have a similar symmetry effect but this band-gap 
is well above Fermi-level and the film energy QSE period 
is given by the free electron result. 

An important remaining issue is whether the bilayer 
periodicity of the freestanding Pb(100) film energy sur- 
vives on adsorption on a semiconductor. A simple es- 
timate of the energy of a thin Pb(100) film grown on 
a semiconductor surface using the "electronic growth" 
model Q indicates that the bilayer periodicity survives. 
In this model, the total energy E t of a metal film grown 
on a semiconductor substrate is given by E t = Eq — E c , 
where Eq is the (relaxed) film energy shown in Fig. 1 and 
E c is the energy gain per surface unit cell due to charge 
transfer at the metal semiconductor interface. This en- 
ergy gain can be estimated in a simple capacitor model: 
E c = C(A$) 2 [l|, |2|, where C is a phenomenological pa- 



4 



rameter specific to the interface and A$ is the difference 
between the work functions of the clean substrate and 
the film. For example, E c can be estimated in the case 
of Si(lll)-(7 x 7) semiconductor substrate by using the 

° 2 

same value for C = 0.033/(eVA ) as used by Yeh et al. 
for similar Pb(lll) films (similar charge density, work 
function, surface energy, etc), and using the measured 
value 4.75 eV for the substrate <I> and our calculated 
values for the film $. Using this estimate for E c and the 
calculated values for Eq, the resulting values for E t are 
0.47, 0.65, 0.54, 0.57, 0.54, 0.59, 0.54 eV for 2, 3, 4, 
5, 6, 7 and 8 layers of the Pb(100) film. Because the os- 
cillations of the work functions are in antiphase relative 
to those of Eq, E c partly compensates the strong oscil- 
lation of Eq but its odd-even layer alternation has not 
been destroyed. 

The above analysis does not necessarily suggest that 
Pb(100) films can grow on Si(lll)-(7 x 7). But it does 
illustrate that the total energy can have an odd-even os- 
cillation when a Pb(100) thin film is grown on a semicon- 
ductor substrate. However, one has to choose a substrate 
with a good lattice match so that the elastic energy is 
negligible. The Pb(100) film should then exhibit "elec- 
tronic growth" if care is taken to relieve possible kinetic 
limitations. 

Finally, we would like to discuss the nature of the quan- 
tization of the QW states as obtained from a scrutiny 
of their wave functions (WF). It is found that the cr- 
and the 7r-like QW states form two independent series of 
states 0. As expected, the lowest QW state derived 
from the bulk s-band has zero nodes, while the second 
one has one node, etc. The bulk s-band accommodates 
20 QW states for a 20 layer film with a maximum of 19 
nodes. The low-lying p z -QW states have a character. 
The lowest p z -Q,W state has 20 nodes and constitutes 
the next state in the a series, followed by the p z -QW 
states having 21, 22, 23, 24 nodes, etc. The higher-lying 
Px,y-QW states with tt character are two-fold degenerate 
and form an independent series of QW states from the a 
series of s- and p z -QW states. The lowest-lying tt state 
(Eg) has quantum number one (zero nodes), the second 
tt level (E u ) has quantum number two (one node), etc. 
The separate quantizations of a and tt states are rather 
general. We have investigated the QW states in linear, 
isolated chains of Na, Mg, Al, and Pb. The QW states 
spectra in these linear chains are quite similar to that of 
the Pb(100) film at T point. The two-fold degenerate tt 
states are higher in energy and quantize independently 
from the a states. 

In summary, we have shown from density functional 
calculations that Pb(100) thin films exhibit quantum size 
effect with an odd-even-layer alternation of film energies, 
film relaxations, etc, which are due to the different sym- 
metry of the stacking geometry in odd and even layer 



Pb(100) thin films. The quantum well states at the film 
Brillouin zone boundary have a two-fold degeneracy in 
even-layer films, which are lifted in odd-layer films. Fur- 
thermore, a large band-gap at the Brillouin zone bound- 
ary in Pb(100) films makes the odd-even alternation of 
the splitting of the bands embracing the Fermi level sig- 
nificant for the film energy. This alternation of the film 
energy is argued to survive on a semiconductor substrate. 
Even-layer Pb(100) films should be especially stable and 
may show up as "magic thickness" . Furthermore, we 
show that the quantum well states in a simple-metal film 
can be classified into a and tt bonded states in the per- 
pendicular direction of the film, which quantize indepen- 
dently. 
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